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ABSTRACT 

Ultra-low emittance (tens of nm) beams can be generated in a plasma accelerator using ionization injection of 
electrons into a wakefield. An all-optical method of beam generation uses two laser pulses of different colors. 
A long-wavelength drive laser pulse (with a large ponderomotive force and small peak electric field) is used to 
excite a large wakefield without fully ionizing a gas, and a short-wavelength injection laser pulse (with a small 
ponderomotive force and large peak electric field), co-propagating and delayed with respect to the pump laser, 
to ionize a fraction of the remaining bound electrons at a trapped wake phase, generating an electron beam 
that is accelerated in the wake. The trapping condition, the ionized electron distribution, and the trapped 
bunch dynamics are discussed. Expressions for the beam transverse emittance, parallel and orthogonal to the 
ionization laser polarization, are presented. An example is shown using a 10-/rm CO 2 laser to drive the wake 
and a frequency-doubled TiiALOs laser for ionization injection. 
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1. INTRODUCTION 

Plasma-based acceleratorJi^ can produce extremely large accelerating gradients, enabling compact sources of 
high-energy beams. Rapid experimental progress has occurred in the held of laser-driven plasma accelerators, 
and electron beams accelerated to multi-GeV energies have been demonstratecPI^ using an intense laser driving a 
plasma wave in cm-scale plasmas. Generation of these electron beams relied on self-injection from the background 
plasma in highly-nonlinear plasma waves.® In the regime where background plasma electrons are self-trapped, 
experiments show that beams with sub-micron normalized transverse emittance can be produced.®® 

To reduce the electron beam emittance further it has been proposed to use a laser to ionize electrons at a 
trapped phase in a plasma wake that is independently excited by a particle beairPor an intense laser.® The laser- 
based methocP^ relies on two laser pulses of different colors: a long wavelength pulse, with large ponderomotive 
force and small peak electric held, excites a plasma wake without fully ionizing a high-Z gas; a short-wavelength 
injection pulse, with small ponderomotive force and large peak electric held, co-propagating and delayed with 
respect to the wake drive laser, ionizes a fraction of the remaining bound electrons at a trapping phase of the 
wake, generating an electron beam. Figure [1] illustrates the basic principle of the two-color ionization injection 
method. This two-color, two-pulse ionization injection concept was hrst proposed in Ref. [9l An additional 
numerical study was performed in Ref. [10] subsequent to Refs. |H1|9| 

In this paper we discuss the trapping condition, the distribution of ionized electrons, and the dynamics of 
the trapped electron beam. We present expressions for the transverse emittance, parallel and orthogonal to the 
laser polarization, of the trapped electron beam, valid for beam or laser wakeheld drivers.^ For the all-optical, 
two-color ionization injection method, it is natural to consider a 10-/im GO 2 laser as the long-wavelength drive 
pulse, generating the plasma wake, followed by a TirALOs laser (frequency-doubled, with 0.4 /rm wavelength) 
for ionization injection. Progress in GO 2 laser technology has opened the possibility of sub-ps pulse durations, 
that would enable efficient (i.e., resonant, with duration of order the plasma period) plasma wakeheld excitation 
at plasma densities ^ 10^® cm“^, and such laser systems are expected to become available in the next several 
years.E^ We present an example of two-color ionization injection using a short-pulse, GO 2 drive laser pulse and 
a frequency-doubled, TiiALOa injection laser pulse. 
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Figure 1. (a) Drive laser vector potential oo, ionization laser vector potential ai, excited plasma wakefield Ez/Eo (red 
dashed curve), and trapping condition 0 < 0 — ())min — 1 + l/7p (blue curve), (b) Drive {El) and ionization {Ei) laser 
electric fields (normalized to the peak of El) and the fraction of the high-Z gas state ionized fi (red dashed curve). 

2. TRAPPING PHYSICS 

2.1 Trapping condition 

The two-color ionization injection concept relies on a large plasma wave (or wakefield) excited by a long- 
wavelength drive laser pulse, with wavelength Ao = 27r/A;o and normalized vector potential amplitude ag = 
eAo/rrieC^, where e and m-e are the electron charge and mass, respectively, and c is the speed of light. For 
efficient wakefield generation the drive pulse should be approximately resonant with the plasma, with kpLg ^ 1, 
where Lq is the drive laser length and ujp = kpC = (dTre^Ue/m-e)^^^ is the plasma frequency, with Ug the electron 

number density. Co-propagating and delayed with respect to the drive pulse is an ionization laser pulse, with 

wavelength A^ = 27r/fcj < Ag and normalized vector potential amplitude ai = eAilmgC?. The ionization pulse 
electric field Ei = (^-nragC?je)ail\i is sufficiently large, Oi/A^ > aojAo, to ionize remaining bound electrons in 
the gas, that may be trapped in the wakefield. 

Trapping requires that the electrons are ionized at the proper wake phase and that the wake amplitude is 
sufficiently large. Behind the drive laser pulse, and assuming af{ipi) <C 1, this condition may be expressed aJ^^ 

1 -7p ^ (1) 

where = kp^ = kp{z — ct) is the phase position in the plasma wakefield of the ionized electron (initially at 
rest), with jp is the Lorentz factor of the plasma wave phase velocity, and (j) is the potential of the wakefield 
(normalized to rrigC^je), with (p^in in the minimum amplitude of the potential. Here we have assumed the 
wakefield near the axis is approximately described by the ID nonlinear wake potential,® whose extrema satisfy, 

cPmin/m..= {Ez/Eo)y2±l3p[[l + {EjEo)^/2] -l| , (2) 

where Ez/Eq is the peak of the accelerating wakefield normalized to Eg = meC^kp/e. The optimal ionization 
phase ipi for trapping is at the peak of the wake potential, <p{il)i) = pmax and Ez{ipi) = 0. At this phase location. 
































the trapping condition is 1 — 7^ ^ < (pmax — In the limit 7p ^ 1, the minimum wakefield amplitude required 
for trapping is 

E,/Eo>{V5-2)^/^:^OA9. (3) 

The longitudinal wakefield amplitude excited by a resonant, circularly-polarized, Gaussian laser pulse is E^/Eq ~ 
(7r/2e)^/^ag(l -I- such that the required laser amplitude for trapping is oq > 0.94. 

Figure [1] illustrates an example of the concept of two-color ionization injection. Figure [TJa) shows the drive 
laser vector potential uq (pulse centered at the origin), the field of the excited plasma wake E^/Eq (red dashed 
curve), the ionization laser vector potential [pulse centered at kp{z — ct) ~ —9.5] Ui <C oq, and the trapping 
condition Eq. o, where 0 < ^ — ^min — 1 + l/7p indicates an electron ionized at rest will be trapped (blue 
curve). Figure jljb) shows the electric field of the drive and ionization lasers (normalized to the peak drive laser 
field), and the fraction of the ionized gas state fi (red dashed curve). The normalized parameters shown in 
Fig. m correspond to the case of a 10-^m wavelength (CO2) drive laser, with oq = 2 (linear polarization) and 
0.47 ps duration (FWHM laser intensity), propagating in a uniform gas, producing a plasma with electron density 
Ue = 10^® cm“^, and resonantly exciting a wakefield. The ionization pulse, with = 0.13 (linear polarization), 
\i = 0.4 ^m (frequency-doubled Ti:Al203 laser), and 120 fs duration (FWHM intensity), ionizes Kr®+ to Kr®+ 
(ionization potential Ui = 230 eV) at a trapped wake phase. 


2.2 Trapped charge 

Tunneling ionization will determine the distribution of electrons ionized, and, hence the trapped charge and 
emittance. The transverse phase space distribution of the ionized electrons was derived in Ref. 1111 The rms 
radius of the transverse distribution of ionized electrons iiP 


where Wi is the ionization laser spot size 
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where re = jmet? is the classical electron radius, a is the fine structure constant, and Ui is the potential of the 
state of the gas used for ionization injection, normalized to Uu ~ 13.6 eV. The parameter A^ is the normalized 
laser field amplitude Ei oc and A^ <C 1 is satisfied at the ionization threshold. (Here we assume the 

Keldysh parameter = [2C/i/(meC^af)] satisfies < Ij such that tunneling ionization is the dominant 
ionization mechanism.) 

For an ionization laser pulse located at the proper phase of a plasma wave of sufficient amplitude E^/Eq, 
given by Eq. o, electrons ionized at rest will be on trapped orbits, and the amount of trapped charge will be 
determined by the number of ionized electrons: Nt = 2'Kax<Jy(-ifing, where fi is the ionization fraction, Ug is the 
number of ions with electrons in the proper charge state, and li is the interaction length. For interaction with a 
single species gas, Ug = n^jZ, where rig is the electron density and Z is the number of electrons ionized by the 
drive laser (or pre-ionized). If the length of the high-Z gas region is longer than the ionization laser Rayleigh 
range, then the total charge will be limited by the diffraction of the ionizing laser, ti ^ Zm = irw'^/Xi, and, 
hence, we expect the charge to scale as Q = eNt ~ e{TTwfA‘^)£ifing oc wf with the interaction length limited by 
diffraction. As shown below, the emittance scales as e oc wf, and, hence, there is a trade-off between increasing 
the trapped charge {Nt oc wf, limited by diffraction) and reducing the emittance. 


3. EMITTANCE 

The achievable beam emittance using ionization injection into a wakefield, driven by a charged particle beam or an 
intense laser, is determined by the initial phase space distribution of the trapped electrons and the plasma focusing 
forces.l^ The rms of the momentum distribution of the ionized electrons, in the plane of laser polarization, is 
approximateljffl^ 


o-p. - a^A, 


( 6 ) 




where A is given by Eq. ([5]). For a small ionization injection laser amplitude of <C 1, the momentum gain 
from the ponderomotive force of the ionizing laser may be neglected. Using Eqs. ([4]) and ([ 6 ]), the initial thermal 
emittance of the beam generated by ionization injection ijA^ 
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Equation © is independent of the driver of the plasma accelerator (e.g., laser or particle beam) and is determined 
from the ionization physics, i.e., the distribution of electrons produced via tunneling ionization. Unless injected 
matched, the emittance of the trapped beam will grow due to phase mixing in the plasma wakefield. 


3.1 Trapped particle dynamics 

An electron trapped in the plasma wakefield, i.e., ionized at the proper phase to satisfy the trapping condition 
Eq. O, will begin to accelerate and rotate in transverse phase space. For an electron with initial transverse 
position xg and normalized momentum ug, the transverse position and momentum will evolve as 

X ~ [a;ocos(fc/3z) + (uo/A:/3o)sin(fc/32)] = cos(fc/3Z + y;), (8) 

Ux = PxIrUeC ~ 7^/^^ [ug COs(fc/3z) - Xgkpg sin(fc/3z)] = -^^^'^kpgrpg S\TL{kpZ + (p), (9) 


where rpg = [xg + [ug/kpg)'^]^/'^ is the initial betatron amplitude, the energy 7 (z) is increasing as the electron is 
accelerated in the wakefield, kp oc 7 “^/^ is the betatron wavenumber in the wakefield, and kpg = kp{z = 0 ) is the 
initial betatron wavenumber. For laser drivers with parameters satisfying ao(l + Og)“^7^ > kpWg/2 (assuming 
circular polarization), or particle beam drivers, the driver will produce a co-moving ion cavity with kp ~ kpjy/^. 

Averaging the betatron amplitude over the initial beam distribution yieldJ^^^ 
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in the plane of ionization laser polarization, and (r^g) = wlIX? 12 orthogonal to the plane of laser polarization. 

If the ionization injection region is short (due to a finite high-Z gas region or short ionization injection pulse 
Rayleigh range), followed by a post acceleration region that is sufficiently long to allow betatron phase mixing, 
then the normalized emittance will saturate to 

e = {{x^){ul) - {xux)'^] = kpg{rlg)/2. (11) 

The initial beam phase-space distribution Eq. m will determine the final saturated emittance Eq. (ED- Matched 
injection in the plane of laser polarization, where the thermal emittance Eq. © equals the saturated emittance 
Eq. (ITT]) , occurs for laser-plasma parameters such that kpgWi = \f2ai. If continuous ionization injection occurs 
over the length of the plasma accelerator, the saturated emittance will asymptote to e = fc/ 3 g(r|g)/\/ 3 . 


4. TWO-COLOR IONIZATION INJECTION EXAMPLE 

In this Section we present an example of beam generation using two-color ionization injection with a short-pulse 
CO 2 laser driving the wakefield and a frequency-doubled TiiALOg ionizing laser. Modeling was done using the 
iNF&RNc(i3 and wariE^PIC codes of blast (Berkeley Lab Accelerator Simulation Toolkit).!^ Figure [5] shows 
the evolution of the transverse normalized emittance of the trapped electron beam versus propagation distance 
z. In Fig. [21 a 10 /rm-wavelength, 10 J, CO 2 laser with a = 2 (linear polarization), 0.47 ps duration (FWHM 
laser intensity), and 155 /im spot size {Zr = 7.5 mm) propagates in Krypton gas, ionizing the gas up to Kr®+, 
producing a plasma with electron density 10^® cm“® (gas density of Ug = 1.25 x 10^® cm“®) and driving a wake. 
A frequency-doubled (0.4 p,m wavelength) TirALOg laser with Ui = 0.14 (linear polarization), 120 fs duration 
(FWHM laser intensity), wt = 20 pm spot size {Zm = 3 mm), delayed with respect to the CO 2 laser, ionizes 
Kr®+ to Kr®+ {Ui = 230 eV) and generates a trapped electron beam. The solid black and red solid curves are the 








Figure 2. Example of two-color ionization injection using CO 2 drive laser: A (10 pm wavelength) 10 J, CO 2 laser with 
a = 2 (linear polarization), 0.47 ps duration (FWHM laser intensity), and 155 pm spot size {Zr = 7.5 mm) propagates 
in Krypton gas, ionizing the gas up to Kr®"*", producing a plasma with electron density 10*^® cm“® and driving a wake. 
A frequency-doubled (0.4 pm wavelength) TiiAEOs laser with = 0.14 (linear polarization), 120 fs duration (FWHM 
laser intensity), Wi = 20 pm spot size {Zm = 3 mm), delayed with respect to the CO 2 laser, ionizes Kr®”*" to Kr®^ 
{Ui = 230 eV) and generates a trapped electron beam. Evolution of the transverse normalized emittance of the trapped 
electron beam versus propagation distance z{mm) is plotted: (black curve) emittance in the laser polarization plane from 
PIC simulation, (dashed black line) e|| theoretical saturated emittance in the laser polarization plane Eq. (I12II . (red curve) 
emittance orthogonal to the laser polarization plane from PIC simulation, and (dashed red line) ex theoretical saturated 
emittance orthogonal to the laser polarization plane Eq. 113}. 


normalized transverse emittance of the trapped electron beam, parallel and orthogonal to the laser polarization 
plane, respectively, calculated using the inf&RNO module of blast. 

For a small ionization injection laser amplitude a? <C 1, propagating in a uniform gas jet, the final normalized 
transverse emittancei^^^ are, using Eqs. (HH), (nni), and ®, 
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where e|| (ex) is the transverse normalized emittance parallel (orthogonal) to the plane of ionization laser polar¬ 
ization. Shown in Fig.[2]are the theoretical predictions: (dashed black line) ey emittance in the laser polarization 
plane Eq. (IT^ and (dashed red line) ex emittance orthogonal to the laser polarization plane Eq. (IT^ . For these 
drive laser-plasma parameters, fc/3o/(^p/’\/2) — 0.9. Good agreement is achieved between the PIC simulation 
and the theoretical estimates Eqs. m and (HU). Phase mixing is not complete after 5 mm, and the emittance 
orthogonal to the plane of laser polarization is still varying. The charge contained in the beam modeled with 
the PIC code is 27 pC. This approximately agrees with the order-of-magnitude estimate of the trapped charge 
over a Rayleigh length, assuming Ri 1, 


Q = eNt ^ engZRi'KWiA'^, 


(14) 


which, for the parameters of Pig. [5J yields Q ~ 34 pC. 

In general, a trade-off exists between the reducing the emittance (e oc wf) and increasing the trapped charge 
{Q oc wf, with interaction length limited by diffraction). Eigure |3] shows the theoretical estimates for the 
beam emittance [parallel, Eq. (fT^ (dashed black curve), and orthogonal, Eq. (IT^ (dashed red curve), to the 
ionization laser polarization] and the charge, Eq. (|14[) (dashed blue curve) versus ionization laser spot size (all 
other parameters are the same as Fig. [^. The points in Fig. [3] are PIC simulation results (calculated using the 
INF&RNO module of blast) for the emittance in the laser polarization plane e||[nm] (black) and orthogonal to 
the laser polarization plane ex[nm] (red), and the charge Q[pC] (blue). 













Figure 3. Emittance [parallel, Eq. (1121) (dashed black curve), and orthogonal, Eq. (11311 (dashed red curve), to the ionization 
laser polarization] and the charge Eq. (1141) (dashed blue curve) versus ionization laser spot size (all other parameters are 
the same as Fig. n. Points are PIC simulation results for the emittance parallel ey [nm] (black) and orthogonal e_L[nm] 
(red) to the laser polarization, and the charge Q[pC] (blue). 

5. SUMMARY AND CONCLUSIONS 

Ultra-low emittance electron beams, on the order of tens of nm, can be generated in a plasma accelerator 
using ionization injection in a wakefield.^l^Eil In this paper we have reviewed the trapping physics, distribution 
of ionized electrons, and the trapped beam dynamics. Expressions for the achievable beam emittance using 
ionization injection into a plasma wakefield were presented.^^^ Note that these results are independent of the 
Wakefield driver (lasei^^ or particle beairPI). An experimentally-relevant example was presented of ultra-low 
emittance beam generation by two-color ionization injection using a short-pulse, CO 2 laser driving the wakefield 
and a frequency-doubled, TitAUOa laser ionizing a Krypton gas. In this example the interaction region was 
limited by the Rayleigh range of the ionizing laser. In general, a trade-off exists between the reducing the 
saturated emittance (e oc wf) and increasing the trapped charge {Nt oc wf, with interaction length limited 
by diffraction). Self-consistent numerical modeling of the beam injection for this example was performed and 
compared to analytic estimates, showing good agreement. It should be noted that this example is not fully 
optimized, and, depending on the laser-plasma parameters, other gases may be considered (e.g., Ne) for improved 
performance. 
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